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We study black hole solutions in R 4 x S 1 space, using an expansion to fourth order in the ratio 
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I. INTRODUCTION 

o : 

If the topology of space-time is i? 4 x S 1 then the only known exact solution representing a black object is the 
fS| . uniform black string with horizon topology S 2 x S 1 Q] . Though this solution exists for all values of the mass it is 
^ ' unstable below a critical value, -Mql, as shown by Gregory and Laflamme 0]. 

Horowitz and Maeda Jj] argued that a uniform black string cannot change its topology into a black hole in finite 
\^ ' affine time, making the possibility of such a transition questionable. They suggested the possibility of a transition to 
jy-j , a nonuniform black string. Gubser jj| showed the existence of nonuniform black string solutions. The non-uniform 

■ black string solution in 6 dimensions was investigated numerically 5] for a range of the mass values above A/ql- 

Nonuniform black string configurations do not exist for masses below the Gregory-Laflamme point, in the region 
CNJ . where the uniform black string solution is unstable [f|. A natural candidate for a black object in this mass range 
is a black hole. Unfortunately, no exact black hole solutions are known in a 5 (or more) dimensional space with a 
compactified dimension. Still, on the basis of physical intuition, a black hole solution should exist for very small values 
of the mass. When the radius of the horizon is much smaller than the size of the compactified dimension, i.e. when 
M — > 0, the black hole should be unaware of the compactification. Then a Myers-Perry |f| solution should become 
asymtotically exact. Indeed, a numerical solution, extending the Myers-Perry solution to larger values of the mass 
was found by Harmark and Obers 0, @ and Kudoh and Wiseman ■ 

Using general arguments Kol ^3 suggested that the black hole branch and the nonuniform black string branch meet 
at a point when the black hole fills the compact dimension. Further arguments for such a transition were presented 

m Boa 

Recently, we have studied a related problem, namely the existence of black holes in Randall-Sundrum [jjj 13 
theories. We used |l5lll6j| an approximation scheme based on the expansion of solutions in the ratio of the radius of 

■ the horizon of the black hole to the ADS curvature. In this paper we employ a similar strategy by expanding the metric 
and other relevant quantities in the ratio of the two natural lengths associated with a black hole configuration, (i) the 

J^j • five dimensional Schwarzschild radius, (A, associated with the mass ^ — ^J&G§M / (37r) and (ii) the compactification 
rS \ length, L, defined as the proper circumference of the compact dimension in the region far away from the mass. The 
dimensionless ratio of these two quantities serves as an excellent expansion parameter. As the solution is an even 
function of \i and L we use e = ( j -) as our expansion parameter. Such an expansion for ADD black holes has recently 
been proposed and Harmark |17l| and Gorbonos and Kol [lg] evaluated the leading term of the expansion. 

To find an unique solution we must fix boundary conditions, both at infinity and at the horizon of the black hole. 
We find solutions in two different regions. In the asymptotic region, r » /i aii 'asymptotic solution', and in the 
near horizon region, r <C L a 'near solution' is found. Using the asymptotic solution we satisfy asymptotic boundary 
conditions but the boundary conditions at the horizon cannot be satisfied. The near solution suffers from the opposite 
problem, as it cannot be used to satisfy asymptotic boundary conditions. To solve those problems we combine the 
two solutions. We match the asymptotic solution with the near solution in the region /i <C r <C L where both are 
valid. This method was used in 1151 to find small black holes in Randall-Sundrum scenario, and in Il3| to find 
small black holes on cylinders. The work in this paper is parallel to and agrees with previous results [rilll£|. but here 
we calculate the metric up to fourth order in /i/L (second order in e). The organization of this paper is as follows: 
In section [H] we describe the general method of our calculations, and list the parameters that will be calculated from 
the metric. In sections 11111 V II we present the detailed calculations up to second order. In section IV 111 we summarize 
our results for the metric, up to second order, and discuss the possible scenarios for black holes of increasing mass. A 
reader not interested in technical details should read the next section, (|Hj) and then proceed directly to the summary 
and discussions (|VII|) . 
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II. THE GENERAL METHOD 



As we indicated in the Introduction we will investigate black hole solutions of the Einstein equation in 5 dimensional 
space when one of the dimensions is compactified. Recently, this problem has attracted the attention of several 
groups. Harmark and Obers[M.ll2j introduced the relative tension of black holes as an order parameter, wrote down a 
generalized Smarr formula [19(, investigated the phase diagram for black objects, and studied the solutions numerically. 
Gorbonos and Kol [Isj proposed an analytic approximation scheme based on the expansion in the ratio of the radius 
of the horizon and the compactification length. This method is similar to the expansion method we used [HEil to 
investigate black holes in the Randall-Sundrum [T3L scenario. We will follow a similar path here and calculate 
further terms of the expansion. 

By necessity, the perturbation method is applied in two overlapping regions. In the asymptotic region, r 3> /i, /i is 
the smallest scale. Therefore, an expansion in \ijL is actually an expansion in fi and in every order the metric is a 
general function of the dimensionless coordinates x/L, 

9ap(x) = ^/Lrg^(x/L). (1) 

n=0 

We refer to this solution as the 'asymptotic solution.' 

In the near horizon region, r <C L, L~ x is the smallest scale. Therefore, an expansion in \ijL is actually an expansion 
in L^ 1 . In every order, the metric is a general function of the dimensionless coordinates x/fx, 

g a p{x) = Y J {n/L) n 9«p%{x/n)- (2) 
We refer to this solution as the 'near solution.' 



A. The equations 

In each region, when we calculate the n-th order terms, Einstein's equations is solved in terms of two functions: 
A gauge function and a wave function, which satisfies a linear, (inhomogeneous), master equation. The differential 
operator for the master equation in the asymptotic region is 

£P_ + 2d_ + 

dr 2 r dr dw 2 ' 

where w is a coordinate in the compact dimension and r is a radial coordinate in 3-dimensional space. In the near 
region the operator is 

^-^- l Hm) + W- iM2 *M-- 3 - (4) 

where R is a radial coordinate and ip is an angular coordinate in 4-dimensional space. The driving term in the master 
equation for the nth order contribution depends on the solution in lower orders. It determines an inhomogeneous 
solution in the wave function. The additional homogeneous solution introduces new parameters that should be fixed 
by boundary conditions. 



B. The boundary conditions. 

In empty space, or for the uniform black string configuration, the cylinder is invariant under translations in the 
compact direction. When a point mass is introduced translation invariance is broken, but a Z2 reflection symmetry 
remains around the location of the mass in the compact direction. So, if we endow the compact direction with the 
coordinate w € [-L/2, L/2] and put the mass at w = then the solution is symmetric about w — > —w and periodic 
in w with period L. 

Asymptotically (far away from the mass), we assume that the metric is Minkowski. The leading order correction to 
the Minkowski metric is of order r _1 . Such a term provides, among others, the four-dimensional Newtonian potential. 

We must impose constraints on the near solution, as well. The horizon of a small black hole has S* 3 topology. We 
require that there is no black string attached to the horizon. Consequently the Kretchmann scalar must be regular 
on the w axis everywhere, when w 7^ 0. Furthermore, the surface gravity must be constant, otherwise the horizon is 
not regular pfj| . 
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C. Matching the asymptotic and near solutions 

To match the asymptotic and the near solutions we must work in the intermediate region p <C p <C L, where 
p is a four-space-dimensional radial coordinate. In this region the functions 9^ y ™ P (^/^) can be expanded in the 
coordinates. Owing to the S 2 x S 1 symmetry the components depend only on p and i/j, the angle between the 
compact direction and the four-dimensional sub-manifold. Then one can use a double expansion in p/p and p/L to 
write the metric as 

oo oo / \ 2n 2k 

ft*Cft0 = ££'& fc, M(£) (£) > (5) 

The expansion includes only even powers due to the ^-symmetry as will be explained later. 

Globally, our expansion parameter is e = (p/L) 2 . In the two regions e has a different meaning. The asymptotic 
solution is constructed as an expansion in p, such that (p j ' p) 2n (p / 'L) 2k = e n (p/ L) 2 ^ n \ so one keeps n constant 
and solves the equations for all values of k. For example, Minkowski metric corresponds to the zeroth order (n = 0, 
k = . . . oo) approximation to our problem, with appropriate coefficients, g^a ■ The first order, n = 1 ,k = 0. . . oo, 
provides linearized gravity (Newtonian potential), etc. 

The near solution is constructed as an expansion in such that (p / p) 2n (p / L) 2k — e k (p/ p) 2 ( n ~ k \ so at fixed k 

we solve for all values of n. Here, k = 0, with appropriately chosen coefficients, g^a , corresponds to the Myers and 
Perry five dimensional black hole (MPS) [6|- The matching is done in the intersecting points denoted by circles, as 
can be seen in Fig.Q). For example, first order of the asymptotic solution intersects zeroth order of the near solution 
at (n — 1, k = 0). We use this information to fix parameters in the asymptotic solution. The calculation is carried 
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Jvlinkowski space 
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'econd order near solution 
Myers-Fe ny First-order ne ar solution 

FIG. 1: (p/ p) 2n (p/ L) 2k terms as points in the (n, k) grid. The horizontal (vertical) lines describe the asymptotic (near) 
solution in increasing order of e. The matching is done at the intersecting points denoted by circles. 

out order by order. Both solutions (asymptotic and near) must be known up to order e" _1 before calculating order 
e n . Within order e n , one should calculate the asymptotic solution first and then the near solution. The sequence of 
calculation of terms in the two solutions is depicted in Fig. (01 • 
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Order of £=u, 2 /L 2 



\x 2 /V 



u, 4 /L 4 



Asympt. soln. 



1, p 2 /L 2 , p 4 /L 4 ,. 



L 2 /p 2 , 1, p 2 /L 2 ,. 



L 4 /p\L 2 /p 2 , 1,. 
■* — -0- 




Near solution 



1, \i Ip , [x /p ,.. 



p 2 /u. 2 , 1, fx 2 /p 2 ,.. 



p 4 /|A 4 , p 2 /|i 2 , 1,. 



FIG. 2: In the intermediate region the asymptotic solution is expanded in p/L, while the near solution is expanded in p/p. 
The table shows the sequence of calculating terms in the asymptotic and near solutions. Terms connected by an arrow are 
identified by the double expansion. A similar table has been presented in [l8| 



D. Physical parameters 



Two parameters can be measured in the asymptotic region, the ADM mass and the relative tension around the 



compact dimension. 8] Expanding the asymptotic metric around the Minkowski metric we write g^y — T]^ 
Then the ADM mass, M, and the tension, b, are defined as 



h.. 



M = 



1 



lim r 2 



4G 5 r 

b = I d i xT ww = 



L/2 



2, 
— n r 

r 



L/2 

1 

4G*5 r^oo 



dw , 



lim r 2 









rL/2 


'2 








dw 


/-L/2 


r 





(6) 
(7) 



It is useful to define the relative tension, n = —b/M. In [8j the pair (M, n) is used to describe a phase diagram of all 
possible black objects. 

The near solution provides the thermodynamic parameters of the horizon, i.e. entropy and temperature. The 
entropy is defined as the three dimensional 'area' of the horizon divided by 4Gs. The temperature is defined through 
the zeroth law of black holes thermodynamics as the surface gravity divided by 2tt. A Smarr formula 0, 0] relates 
the entropy, temperature, mass, and relative tension as Q 



TS = 



2 - n{M) 



M 



(8) 



III. ZEROTH ORDER IN e 



In the asymptotic region zeroth order in e means that we neglect the mass. Therefore, the configuration is asymp- 
totically flat (since the mass is confined to a small region). So, the asymptotic form of the metric is Minkowski on a 
cylinder 

ds 2 = -dt 2 + dr 2 + dw 2 + r 2 (d0 2 + sin 2 6W0 2 ) , (9) 

where w is the compact coordinate such that w € [— L/2, L/2). 

In the near region zeroth order in e means that L — ■> oo. The mass point is in five dimensional asymptotic flat 
infinite manifold. Therefore, the near solution is given by the Myers-Perry Q solution (MPS), 

ds 2 = -(1 - 4)* 2 + — ~~~2" + P 2 d^ 2 + p 2 sin 2 ^(d6 2 + sin 2 6d<t> 2 ) , (10) 
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where p is the five dimensional Schwarzschild radius. In order to match the metrics (0 I10[l one needs to specify 
a coordinate transformation between (t,r,w,6,(f>) and (t, p,ip,6,(f>). We choose to work with a specific coordinate 
transformation throughout the calculation of all orders. (One can, however, choose to redefine the transformation in 
each order). The coordinates t, 9, and <fi are shared by the coordinate systems. Furthermore, we want to preserve the 
periodic structure in w so we choose the transformation 

L . ttw , 

pcosip , iu £ |— L/7T, L/7rJ , (Ha) 

(lib) 

7TT* 

(lid) 



u = 


— sin — — 

7T L 


r = 


psaitp , 


P = 


•\/r 2 + v 2 


^ = 


r 


arctan — = 



Lshx(mv/L) 

One can verify that the metrics (|9I10|I agree up to order /jP(1/L)°, 

We show now that the expansion contains only terms with integer powers of (p/L) 2 = e. It is sufficient to establish 
that the solution must be an even function of p and of L (at fixed p 2 ~ MG5). The solution is an even function of p, 
because the driving term, the zeroth order near metric (|10f> depends only on p 2 . The zeroth order near metric, l|l(Jf) . 
combined with I jlljl is also an even function of L. The conditions imposed on the metric also force it to be even in 
L. The mass is located at p = 0, (r = w = 0). The configuration of a mass point on a cylinder has a Z2 symmetry 
about w = 0, {ip — 7r/2), requiring that the terms of the metric depend on w and L through functions of the form 
cos(2nnw / L) . 



IV. FIRST ORDER IN e 
A. The asymptotic metric 



In the asymptotic region, r»/i, we expand the metric as 

g a p = WaP + th a p , (12) 

where r] a p is the flat metric (J3J. 

The general form of the solution for h a p appears in many places in the literature . We re-derive it here using the 
matching method. We expand the Einstein equation to first order in e and get a set of homogeneous linear equations 
for h a p. Using the coordinate transformation w — > w + eW(r, w) we choose the gauge h rw = W <r (r, w). Next we solve 
G r w = G ww = for h rr , and G rr = for h ww , where G M „ is the Einstein tensor. The solution is given by 

01 7* 

h r r = n^tt.r ? (13a) 

r 2 

OL T 

h ww = -htt,r + 2Ww(r,w) , (13b) 

r 2 

h rw = W, r (r,w) , (13c) 
hee = h^ = , (13d) 

where the gauge function satisfies the symmetry properties W(r, w + L) = W{r, w) = — W{r, —w) and a is a constant. 
The function hu satisfies the equations 

2 

htt.rr H htt, r + h u ww = , 

r 

h tt (r,w) = h t t{r,-w) ; h tt (r,w + L) = h tt (r,w) . (14) 
The solution to Ea. (|14f> can be written as a Fourier series 

u i \ X 2, exp(-27rnr/L) 

htt{r,w) = > c„ cos{2imw / L) . (15) 

n=o r 

The constants c„ are determined by matching i|15|) to the zeroth order near solution (|1(J|) . Then it follows that up to 
leading order in L _1 the two solutions must be identical, th u = [J? p~ 2 ■ The constants c„ are given by the Fourier 
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expansion of fx p 



2„-2 



c 1 



L/2 



-dw 



e— exp(— 2nnr/L) = — [ 
r L J _ 



r L J_ L /2 r 2 + (L/7r) 2 sin (irw/L) 

2 

cos(2Trnw / L)dw 



c = ttL + O{L 



(16a) 



2 

I -L/2 r 2 + (L/7r) 2 sin 2 (7rw/L) 
Using l|16[) we are able to sum Ea. (|15|l as 



c n = 2itL + O(L ) . (16b) 



ttL sinh^ 

h tt (r,w)- 



r cosh^-cos 2 ^ 



For r>L the metric is given in terms of the 'zero mode,' independent of u>, and can be summarized as 



flu 

Qrr 



-l 



5« 



= l 



2a - ttL 
' 2r ' 
2a + ttL 



2r 



+ 2W, w (r, w) 



eW r (r, iu) 



The constant a, which appears in the metric l|18[l is related to the tension in the compact dimension 



b = / d 4 xT ww = 



4g: 



lim r 



,L/2 


'2 






/-L/2 


r 



dw 



aeL 
2G~ 5 



(17) 

(18a) 

(18b) 

(18c) 
(18d) 

(19) 



The tension b should be zero because we deal here with linearized gravity, with no interaction between the mass and 
its periodic images. Consequently, we must set a = 0. 

Dimensional analysis also requires a = 0. First, we can deduce from Eas. (|13fl that a cx L. The asymptotic solution 
is an expansion in the mass, so a must be independent of the mass and proportional to L, the only quantity of the 
correct dimension at our disposal. As we have discussed earlier, function (|17|l is even in L. If a = cL ^ then the 
functions h rr and ft„, in (|13|) . would have a mixed symmetry under L —* —L. Note that this symmetry argument 
can and will be used in higher (nth) orders of the expansion to eliminate terms of the form a n /r from g rr and g ww . 

Next, we turn to evaluate the conserved ADM mass using Eas. (|18[) 



M = 



e 

4g; 



lim 



L/2 
-L/2 



dw 



3vr/i 2 
^G7 



(20) 



The effective four-dimensional Newton constant is defined by gtt = — 1 + 2G^Mjr. Using Eas. ljlSI 120(1 we find that 
G4 = 4Gs/(3L). We should emphasize that the physically measurable quantities, the relative tension, the mass, and 
Newton's constant, are independent of the gauge function W(r, w). 

As we mentioned earlier, in higher orders of the expansion the relative tension in the compact dimension is not zero. 
However, we require that the ADM mass is completely determined by the first order, such that higher orders do not 
change Ea. (|20|l . This way, we make sure that the expansion in e is also an expansion in M and that thermodynamic 
properties of the black hole are well defined in each order. We will return to this issue in section IVll when we discuss 
the second order contributions. 

The gauge function W(r, w), which appears in (jlHjl . can be partially determined by matching to the zeroth order 
of the near solution. However, we prefer to determine the function completely by matching to the first order of the 
near solution, as well. 



B. The near solution 



In the near region it is convenient to use the coordinates (t, p, tp, 9, <fr) just like in Eq. (fTU|l . We use the following 
ansatz for the metric: 

ds 2 = -B(p, 4>)dt 2 + tfj dp 2 + 2V(p, ib)dpdib + p 2 U(p, i>)d^ 2 + p 2 sin 2 ^{dO 2 + sin 2 d(/> 2 ) . (21) 
B(p,ip) 
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The functions B, A, U, V are expanded to first order in e as 

B(p,ip) = l-^-eBxM) , (22a) 

AM) = 1 + eAxM) , (22b) 

U(p,i/>) = l + eUiM) > (22c) 

V(p,ip) = zVx{p,i>) • (22d) 

The first step is to use a gauge transformation of the form p — > p(l — ei*i), ^ — ► ^ + e-Fi tan-0, which leaves the last 
term of l|2I|) unchanged to leading order, to choose a gauge where 



This gauge choice simplifies the equations for the rest of the functions to be determined. 

Next, one can solve the equation G pl p = for B\. The rest of the equations can be solved in terms of a single 
function, Hi(p,ip), obeying a second order partial differential equation. The solution is 

Bi = , (24a) 

P 



Ai = 4H ^f p ^] + 2F 1 (p, i>) + 2pF 1 , p (p, ff>) , (24b) 
sm ip cos ip 

= _ 2 6 P H lM ) - ( y - f) Hl AP,*) _ 2 , . F| , , , „ 2 !;| „ , , F „ } ,, J( , 
p sm ip cos^ ip 

Vl = , 2 P(2 2 f~2 M ? -[3^i(p»-tan^i^(p»] 

(p — /x z )sin ip cosip 



_ 2p HyM) + ^_^ ^ 
sm cos ^ P ~ A 1 

The wave function TJi (p, ^) satisfies the differential equation 

{p 2 - /i 2 ) (hi, pp - ^#i,p) + H hW - 4 cot 2 ^i,v - 3ffi = . (25) 
The solution of Eq. I|25|) can be written as 

H 1 (p = Rp,ip) = 2Ry/R 2 - 1 J dv[a{v)Pl{2R 2 -l) + b{v)Ql{2R 2 -1)}K{tP) , (26) 

h v {ip) = -^=sin 3/2 2V> [7rcosa^P^ 3/2 (cos2V') - 2 sin o^Q 3 / 2 (cos 2-0) 
V 2vr L 

= (2f - f ) cos[(2i/ + 3)V> - a„] - (2f + 3) cos[(2z/ - 1)V> - a v ] , (27) 

where P^ 1 and Q£J are associated Legendre functions of the first and second kind. The functions a{y) and b(v) 
will be fixed below using the constraints on function H\. These constraints come from symmetries, from regularity 
requirements, and from boundary conditions. 

1. Conditions at ip = n/2 

Z2 symmetry about the direction ip — tt/2 (w = 0) implies the conditions 

ffi,^(p,7r/2) = , (28a) 
Fi,^(p,7r/2) = , (28b) 

Fx (p, tt/2) = -6H 1 (p,7r/2)+ 2/J ~ M H lt( ,(j>,ir/2) ■ (28c) 
It can be verified that condition (|28a|l is automatically satisfied by (|26|) . 
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2. Conditions at ip — 

For a black hole configuration, unlike for a black string configuration, the mass is localized at the origin. Then the 
components of the metric, (|22[l . are finite at ip = and arbitrary p. This will only be true if 

#i0,0) = ffi )V ,(p,0) - ffi,^(p,0) - (29) 

It can be verified that (|29|l implies cos a v = 0. Then l|27|l simplifies to 

h v (ip) = (2v - 1) sm{2v + 3)tp - {2u + 3) sin(2i/ - l)ip , (30) 

where we have omit a possible overall sign. 



3. Evenness in p 



As we mentioned earlier, the expansion in e = (p/L) 2 and the fact that the zeroth order of the near solution 11 0|l 
depends only on p 2 imply that the metric is even in p. This means that the function Hi should be even in R = pj p. 
For large R the Legendre functions in (|26Jl behave as 

RVR 2 -lP,3(2i? 2 - 1) ~ R 2u+2 x f P (R- 2 ) , (31) 

rVR 2 - lQl(2R 2 - 1) - R- 2v x f Q (R- 2 ) , (32) 

where fp/Q are analytic functions. Since we require that H\ is even in R, v should be integer. The integral in l|26|) is 
reduced to a sum over integer values of v. 

oo 

Hi{p = RpL,ip) = [a (R 2 - 1) + b ] sin 3 vb + 2R^ R 2 - 1 ^ [a„P„ 1 (2i? 2 - 1) + b n Q\{2R 2 - 1)] 

n=l 

x [(2n- l)sin(2n + 3)^ - (2n + 3) sin(2n - i)ip] . (33) 

The case n = requires special attention, since the Pq = 0. 



4- Requirement of a Killing horizon 



The metric (1211) is static. Therefore, the surface B — is a Killing horizon. The normal vector i? Q should be null 
on the horizon. This implies that on the horizon = 0, thus the horizon is located at constant p — pn- Using 
metric (|22ll we expand pn in e as pu = fx(l + e(i)- The conditions B(pH,tp) = and B t1 p(pH,ip) = restrict the 
gauge function Fi 

F 1 (p = n,il>) = -C 1) Fi,^(p = A*,V) = 0. (34) 
The surface gravity for metric (|21|l is defined as 



2^4GM0 



(35) 



The surface gravity should be constant on the horizon, otherwise the horizon is singular [2£j ■ When we expand the 
metric in e, the surface gravity must be constant in every order of the expansion. We use metric (|22|) and conditions 
(|3"P to evaluate the surface gravity in order e 

n=l-2e— -f- - = l + exi- (36) 

sin cos ip 

The requirement of a constant surface gravity constrains the function Hi . At p = p the Legendre functions in Eq. Ij33(l 
take the values 2Ry/R 2 - \P^(2R 2 - l)\ R=1 = 0, 2Ry/R 2 - IQ\{2R 2 - l)\ R=1 = -1. We use the representation © 
to evaluate the surface gravity i|3t)|) 

oo 

.\^h-(9.n- -iW>.n + $\ — 

sin ip 



6b - 4 Vfo„(2n - l)(2n + 3) Sm '- 2n + 1 ^ = Xi = constant . (37) 

— ' sm in 



n=l 
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The set {sin(2n + 1)0} is complete on the interval [0, 7r/2], therefore, the solution to Eg. ((37(1 is 

-66 = Xi, K>o = . 



(38) 



In other words, the sum in Eg. 1(33(1 contains Legendre functions of the first kind only. The functions 
2RVR 2 - lPn(2R 2 - 1) are polynomials of order n + 1 in (2i? 2 - 1). 

C. Matching the near and the asymptotic solutions 

The remaining free parameters of the near solution, {a n , bo, Ci}j must be determined from matching the two gauge 

functions, W(r,w) in the asymptotic solution 1(12(1 and F 1 (p,tp) in the near solution 1(221) . We use to transform 

the asymptotic solution to (p, 0) coordinates and then we expand it in L~ l . The gauge function W(r, w) is also 
expanded as 



00 n 

7TW x — ^ / 7T \ zn 



(39) 



n=0 



We keep the explicit factor cos ^ in ((39(1 to insure that W(r, L/2) = 0. 



1. Matching in zeroth order 



We start with zeroth order in L . We transform the asymptotic metric using 111(1 and expand it to zeroth order 
in L~ x . 

,,2 

(40a) 
(40b) 



B A 


= 1 - 


M 

" p 2 ' 


A A 


= 1 - 


fJL 2 cos^(cos^ — 2p 2 Wo, p (p, ifi)) 


9 1 

p 2 


U A 


= u 


p 2 sin 0(sin ip - 2pW ,^ (p, 0)) 


9 ' 

p 2 


V A 


= P 2 


[cos ipW ^(p,ip) - psimpW , p (p,ip)] 



(40c) 
(40d) 

where the superscript A stands for 'Asymptotic' and the functions are defined in ansatz (|21ll . Matching 140(1 to zeroth 
order of the near solution ((10(1 determines the function Wo 



Wo(p^) 



cos0 
^2p~ 



(41) 



2. Matching in 0(L~ 2 ) 



The asymptotic solution is 
1 



B A 
A A 



p 2 p 2 n 2 



p 2 3L 
.2 



= 1 + ^7 



7T 

L 2 

2 



2 (1 - cos 4 0) , 

p 2 (41 cos 2-0 + 8 cos 40 — cos 6-0) 



p 2 cos 4 + 



9G 



+ 2p 2 cos0W 2 , p (p,0) 



U A = 1 + ^ 



7T 
L 2 



/? cos sin 



2 /z 2 sin cos 2 0(5 sin 30 — 37 sin 0) 2 sin0 
, • _ 2/f 
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-W / 2,v.(p,V') 



(42a) 
(42b) 
(42c) 



V A = - 



7T 2 

-2 [p 3 sin cos 3 + p 2 psimpcos 3 0(sin 2 + 2) + p 2 psm%pW2, p (p, ip) — p 2 cosipW2^(p, 0)] . (42d) 

The near solution 1(22(1 should be expanded to second order in p, which appears in e = (p/L) 2 and in the rescaled 
coordinate R = p/p. The large R behavior of the terms of Hi, ((33(1 . is 



2Ry / R 2 - IP^{2R 2 - I) 



2{2n)\ 
n!(n- 1)! 



R 



,2n+2 



(43) 



10 



These expressions contribute by terms of order eHi ~ L~ 2 /i~ 2 ™ in metric (|21|l . We deduce that n > terms contribute 
by negative orders in p and must be eliminated. Therefore we impose a n = if n > 0. For similar reasons, gauge 
function F\ should also be a polynomial of order 2 in p. Thus, functions H\ and F\ are 

ffi - [a (^ 2 p 2 -l) + Msin 3 V , (44) 
Fi = (m~V - - Ci , (45) 

where we have already imposed F\{p — p,tp) = I).i4|l. The exact form of the periodic function /2(V0 will t> e fixed 
below. Then the near metric, expanded to second order in //, is 

= + (46a) 

^ = 1 + ^ [l2aop 2 + 6p 2 /2W-2^ 2 (6ao-66o + / 2 W + C2)] , (46b) 



w = 1 ^ [p 2 [2a + / 2 (^) + cot^W] - M 2 [4ao - 65 + / 2 (^) + C2 + cot^ty)]] , (46c) 



V N = ^^[-4ao-2/ 2 W+cot^W] • (46d) 
Comparing 1)42(1 and (|46[1 we find that 

/aW) = -^(3 + cos2^)sin 2 ^ , (47) 

^ 2 , ^ 

1^2 (p» = 7 h76 + 11527r" 2 6 "57cos2i/'-12cos4i/' + cos6?/'- 1927r" 2 Cil . (49) 

192 cos tp 

At this point the near solution still contains two free parameters, bo and £i- The /^-symmetry condition, i|28c[) . 
imposes one constraint on these parameters 

TT 2 + 6Cl f . 

bo = -gg- ■ (50) 

A constraint, determining £1, is derived from the first law of black hole thermodynamics. It will be discussed in the 
next section. 



D. Black hole thermodynamics 

The zeroth law of black hole thermodynamics states that the temperature of a black hole is 

T= 2^' (51) 

where k is the surface gravity, which is constant on the horizon. To calculate the temperature for the near solution 
we use Eqs.lgSlEniEHI to find that 

r - i ^ e - < 52 > 

The first law of black hole thermodynamics is dM — TdS, where the entropy is proportional to the area of the 
horizon 

Ah i 



S = AG 5 = AG 5 Jo # io de 'o # ^ 9 ^ me9 ^\ P =P H ■ ( 53 ) 



We calculate the entropy for the near solution and find that 
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According to the first law the temperature can be found as T 1 = dS/dM . The mass appears in the entropy only 

(55) 



through (j, = y/8G a M/(3n) and e = (p/L) 2 , so if we combine the first law and the zeroth law we get 

T- 1 = 271741 + 66 e] = 2vr/i[l + 5&e] . 
This fixes the last parameter. If we combine Eas. (|5()l55(l we find that 

5tt 2 



7T 

24 : 



b = 



144 



To summarize, the near metric, to first order in e, is 



B N = 

A N = 
U N = 

v N = 



i-!L- 

p 2 



Yip 2 



[4(l-cos 4 V)(p 2 -^ 2 ) + /i 2 



1) cos 4 ip - 1 



l + ^ 2 e(^ 



1) sin 2 ip cos 2 tp 



2 P s 
-7r e— 5 sin ^1 cos V 



The location of the horizon is at pn = + 7r 2 e/24). The entropy and the temperature are 



S = 
T = 



TT 2 p 3 

2G 5 
1 

2ixp 



1 - 



57T 2 /i 2 

HI 2 " 



(56) 

(57a) 
(57b) 
(57c) 
(57d) 

(58) 
(59) 



These expressions are in agreement with previous results |17l Il8j . The only freedom left in the first order metric is 
the n > 1 terms of the gauge function in the asymptotic solution (|39|l . The n = and n = 2 terms are fixed in the 
region p <C L by our matching procedure, as given in Eas. l41l49|) . In addition, in the asymptotic region, r»i, the 
metric should be Minkowski, therefore, we require that W(r ^> L,w) = ©(r 1 ). A form of the gauge function which 
is consistent with these conditions appears in appendix lAl 



V. HIGHER ORDER CORRECTIONS 



Prior to completing calculations in second order we describe our general procedure for calculating higher order 
contributions. First we consider the asymptotic solution, which is fully determined (up to gauge freedom) by the 
lower order contributions. 



A. The nth order asymptotic solution 

In the asymptotic region we expand the metric in e as 

00 

ffM^V' + E^M • (60) 

71=1 

Assume that we know the solution up to order (n — 1) and we intend to obtain the nth order solution. Einstein's 
equation is linear in the h^ v but due to lower order contributions it is inhomogeneous. The solution is similar to 
(|13f) but, in addition to the solution of the homogeneous equation includes extra terms corresponding to a particular 
solution (the particular solution is denoted by h^), as follows 

h<${r,w) + — -^h%l , (61a) 
r 2 

h%l{r,w)-^-\h ( £l + 2W%Xr,w) , (61b) 

hi n J(r,w) + Wp\r,w) , (61c) 
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where a is a constant, and W^ n \r, w) is a gauge function which is periodic and antisymmetric in w. The function 
h[t is also periodic and symmetric in w. It satisfies the equation 

/>) + £?>) + - (621 

"-ti.rr ' r n tt,r + IL tt,ww ~ ° I l uz J 

where S^™' is the deriving term, which depends on lower orders, /i^), k < n. The homogeneous solution to Ea. i|62|) 
can be written as a Fourier series. Consequently, we have 

h%\r,w) = h ( t f(r 7 w) + f2c k eM ~ 27Tkr/L) cos(2irkw/L) . (63) 

k=Q T 

The particular solution h^ v is completely determined by the lower orders of the asymptotic solution (without any use 
of the near solution). The free parameters are a and the set {cfc}. These can be determined by matching to the lower 
orders of the near solution as follows. Take gu from the near solution up to order g™ -1 and expand it in /i. Take the 
term of order p? n and find its Fourier series just like in Ea. (|16fl . Compare that Fourier series with the Forier series of 
/ijj to lowest order in L _1 , and determine the constants {c^}. 

Just like in first order, we can show that a n = 0. Solution (|61|l should be even in L. The functions hS n ^ and 
are even functions of L since these functions are determined by lower orders but for dimensional reasons a n ~ L is 
odd. Therefore, a n must vanish. 

Next, we have to make sure that the definition of the ADM mass does not change. So, we require that 



lim r 2 



L/2 



-L/2 



iv rr "'WW, 



dw = Q . (64) 



In general, the tension in the compact dimension, (|19f) . does not vanish. As a result, the effective four-dimensional 
Newton's constant, G4, acquires a correction of order e™. This means that G4 depends on the mass and the equivalence 
principle is violated. 

At this point the asymptotic solution is determined up to order e™, except for the gauge function W^ n \r,w), which 
is (partially) determined by matching to the near solution to orders up to e™. 

B. The nth order near solution 

Near the horizon it is convenient to use the metric in form 121|) and expand it in e as follows 

BM) = l-^-Y^e n B n {p^) , (65a) 
AM) = l + 5>"A„(p,0) , (65b) 

n=l 

UM) = l + ^e"LT n ( p ,^) , (65c) 

n=l 

V(p,il>) = ]Te"K(p» • (65d) 

n=l 

We need the near solution up to order (n — 1) and the asymptotic solution up to order n to calculate the nth order 
near solution. Again, the nth order Einstein's equation is a linear inhomogeneous equation in g' n ) = {B n: A n , U n , V n }. 
The inhomogeneity depends on g( k<n \ We solve the equations, in a way, similar to that for the first order correction. 
The first step is to apply a gauge transformation of the form p — > p{\ — e n F n ), ip — > ip + e n F n t&nip to arrive at a 
gauge, in which 



*o(2p 2 - M 2 ) A , 1 \ 1 cot ip 5 
4(p 2 - n 2 ) 



V n (p, i>) = -7775 —A n ^(p, ip) H U njP (p, ip) . (66) 



Next, one can solve the inhomogeneous equation G p ^, = for B n . The rest of the equations can be solved in terms of 
a single function, H n (p, tp), obeying a second order partial differential equation, which is similar to 125|) . 



(p 2 - p 2 ) ( H n . pp - -H n;P I + H n ^, - 4 cot 2ipH n ^ ~ 3H n = S n , (67) 
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where S n depends on the lower order corrections. The solution of Eq. I|67ll can be found by the method of separation 
of variables, just like it has been done when we have calculated the first order correction (|26|) . The eigenfunctions 
of Eq.ljHZJ) are the functions /i„(0), which appear in Ea. (|2*7jl . The driving term, S n , can also be expanded in the set 
h v (ip)- In addition, the boundary conditions at tj) = 7r/2 and at = in nth order are the same as in first order, 
Eas. l28al I29fl . And just like in first order the function H n should be even in p. So, the general solution of Ea. (|67ll is 

H n (R = p/n,iP) = [a (i? 2 -l) + 6 +Po(i?)]sin 3 V 

OO 

+ [ZRVR 2 - 1 (a n P^{2R 2 - 1) + b n Ql{2R 2 - 1)) + p n (R) 

n=l 

x [(2n - 1) sin(2n + 3)<0 - (2n + 3) sin(2n - 1) 0] , (68) 

where {p n (R)} is a particular solution of the inhomogeneous equation (|67[) . The particular solutions are completely 
determined by the lower orders of the near solution (without using the asymptotic solution). The homogeneous part 
should be completely fixed by matching to the asymptotic solution. The matching is done as follows: Take the 
asymptotic solution up to order e n and transform it to the (p, ip) coordinates, using l|ll|) . Expand the solution in L , 
take the term of order L~ 2n and compare it to the near solution of order e n to fix the parameters {a n , b n }, and the 
gauge functions F n (p,ip) and (r, w). At this point the function W^ n '(r, w) is determined only up to order L~ 2n 
and one should fix it completely before proceeding to the next order. In the next section we apply this method to 
calculate the second order contributions. 



VI. SECOND ORDER IN e 



A. The asymptotic solution - second order 

Up to first order, the asymptotic solution is given by equations ©, (|T2l . (TT*5|) . and l(T7|) . The second order solution 
is described in the Appendix lAl It contains an inhomogeneous part, hP>(r, w), which is completely determined by 
the first order solution, and a homogeneous part, which includes the wave function h tt (r,w), and the gauge function 



A2) 
Hi 



1 {H) 

Ht 



l {2) 



.2) 



/,(2) - -LhW +h (2 

n rr 'Ht,r ~ "rr > 

h {2) - --h {2) + 2W {2) + U 2) 



h (2) = W (2) + j l (2) 



(69a) 
(69b) 

(69c) 
(69d) 



The wave function h\f^ satisfies the homogeneous part of Ea. (|62[) . and can be written in the form of Ea. (|63|l 



Ck- 



exp(— 2irkr/L) 



cos{2-Kkw/L) 



(70) 



fe=0 



To fix the constants, {c^}, we follow the prescription that appears after Ea. (|63|) . We take g tt from the first order near 
solution, l|57a|) . and expand it in p, up to order p , to get 



-l + p z 



7T 2 (1 — COS 4 VO 



3L 2 



4 tt 2 (4cos 4 V'- 3) 
1 12LV 



(71) 



Then we take gtt from the second order asymptotic solution and transform it to the (p,i/)) coordinates, using i|ll[l . 
After expanding it in L , up to order L~ 2 , we obtain 



5t1 = -l + M 2 



1 7T 2 (1-COS 4 0) 
~2 ' 



9 



3L 2 



2 cos 2^ + 1 TT 2 (3 cos 60 + 124 cos 4 + 21 cos 20 + 4) 



Ap 4 



96LV 



(72) 
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We find h^P from gfi. — g$ and calculate its Fourier coefficients 



7T 3 £ 

' 12r 



l-2^(6fc 2 -l)e- 27rfcr/i cos 



2irkw 



fe=i 



tt 3 L sinh ^ (*- cosh + 16 cosh cos ^ + 5 cos ^ - 20) 



24r (cosh 2gt - cos 2™)' 



The asymptotic (r 3> L) form of the metric is 



,(2) 



"'WW 



TT 3 L 

12r ' 

7T 3 L 

12r ' 

— + 

Ylr - w 



The contribution to the ADM mass vanishes indeed, 



AM oc lim r 2 



L/2 
-L/2 



-/l (2) - /l (2) 
'"rr ww.r 



TT 3 L 2 

dw = — 2r 2 W {2) 



w=L/2 
w=-L/2 



= 



where the last equality follows from boundary condition l|A6|l . The tension l|19|) is now 



b = / d"xT ww = 



— — lim r 2 

4G 5 r^oo _ 



,L/2 


'2 








dw 


/ -L/2 


r 





(73) 

(74a) 

(74b) 

(74c) 
(74d) 

(75) 
(76) 



B. Second order correction to the near solution 

Following section IVBI we find that the second order solution is 
B 2 ■■ 



%F^) + [3(P 4 - /*V + ^) cos 4^ 



p^ 216/i 2 p 2 
4(4p 4 - 6p 2 p 2 + 3// 4 ) cos 20 - 5(7p 4 - 21 p 2 p 2 - 3/i 4 ) 



A 2 
U 2 

v 2 



sm V cos V' 864/i 4 



[^ 2 (105/9 2 + ll/i 2 )cos6V> 



(77a) 



= -2 



(-64p 4 + 138p 2 p 2 + 58/i 4 ) cos 4^ + (1472p 4 + 1383/i 2 p 2 + 237//) cos 2^ - 768p 4 + 774p 2 p 2 - 2p 4 ] (,77b) 
6pH 2 (p,^) - (2p 2 - p 2 )H 2tP (p^) 



p sin ip cos 2 ip 
p(2p 2 - p 2 ) 



2 tan 2 0F 2 (p,i/>)-2 tan ^F 2 ^ (p, ip) 



(77c) 



(p 2 — p 2 ) sin ip cos ?/> 

^,3 



[3H 2 (fnp) - tanipH 2 ^(pi/>)] 



2p 2 H 2 JpiP) 



sin 2 tp cos ?/> 
TT 4 p(2p 2 - ^ 2 )sin2^ r 



/z 2 (15p 2 + ,u 2 ) cos 6*0 



+ (-8p 4 + 29/i 2 p 2 + 4p 4 ) cos40 + (184p 4 + 185p 2 p 2 + 11/i 4 ) cos2i/; - 96p 4 + 99/i 2 p 2 - 16p 4 ] 
The function H 2 satisfies the differential equation l|67|l with 



(77d) 



S 2 



5tt 4 (9p 2 + p 2 ) , 7r 4 (20p 4 - 71p 2 p 2 - 23/t 4 ) / * \ 7r 4 (308p 4 + 310 M 2 p 2 + 71 M 4 ) 

2 MV>) ^ ' ' 



48384/1 
7r 4 (432 j o 2 + 41p 2 ) 



34560p 4 
23^W- M 2 ) L 

1TO o 4 h 0W 



24192/t 4 



to WO 



17280p 2 yYJ 1728/t 4 

where the functions h n are given in Ea. (|30|l . 



(78) 
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The solution of Ea. H67|) with inhomogeneity l|78[) can be written in form Eq.JBSJ) with p n >4(p) — 0. For each of the 
non vanishing p n (p) one can add combinations of the homogeneous solution which includes the Legendre functions as 
appear in Ea. 1)68(1 . These should be chosen such that p n (p) is regular at p — 1, and is a polynomial in p of the lowest 
order possible (in principle, the homogeneous solution is of order 2n + 2). In fact, the particular solution for Ea. H67|) 
with inhomogeneity 178|) includes polynomials of order 4 in p, only, 



H£(p,i>) = 


Po(p = 


Rp) 


Pi(p = 


Rp) 


P2(P = 


Rp) 


P3(P = 


Rp) 


Pa(p = 


Rp) 



4 
n=0 

23R 4 



Pn(p)h n (lp) , 



13824 ' 
432.R 2 + 41 

138240 ' 
462i? 4 + 156E 2 + 71 
580608 
24i? 4 - 75R 2 - 23 
1658880 ' 
9R 2 + 1 



774144 



(79a) 

(79b) 
(79c) 
(79d) 
(79e) 
(79f) 



In addition, to avoid negative powers of p, the homogeneous part of H2 and the gauge function F 2 should be 
polynomials of order 4 in p = Rp. 



H 2 (p,ip) = H[{p 7 ij)+ [a (i? 2 -1) +6 ] sin 3 V + 8 ai i? 2 (i? 2 -l)(l + cos 2 V) sin 3 V , 
F 2 (p, i>) = (R 2 - l) 2 / 4 W + (R 2 - 1)/ 2 (V) + foW ■ 
We now turn to imposing boundary conditions. 



(80) 
(81) 



1. Conditions at the horizon 



The horizon is located at a constant p = pa- We set ph = 1 + 7r 2 e/24 + e 2 ^ and solve B(pn, tp) = for fo(tp). We 
find 



cos 8^ cos 6^ cos Atp 7 cos 2ip 583 



2304 432 



36 



72 



6912 



(82) 



The surface gravity is constant (due to the fact that the Legendre functions of the second kind are not included in 



K= 1- 



57r 2 e 
24 



113tt 4 
3456 



(83) 



2. Conditions at 'infinity' 

We match the near solution, l|77(l . to the asymptotic solution, l|A8(l . expanded to order L~ A . Comparing the 
expressions for g u , we find that 

7T 4 

fM) = [29 cos 6V> + 54 cos 4^ - 237 cos 2ip - 166] , (84) 

8640 

7T 4 

J 2 {iP) = — — — — — [15 cos 8tp + 188 cos 6ip + 1248 cos Aip + 3156 cos 2ip + 2993] . (85) 
17280 

Matching the other components, as well, we find the constants 

3117T 4 

ao = -13824 ' (86) 
897r4 

ai = 23040 ■ (87) 
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The gauge function, W^ 2 \ is 

= - 



3i 4 cos ip 



n 2 L 2 



- (327r - 64^ + 15 sin 2t[) + 7 sin 4-0 + 5 sin 6ip) 



32p 3 1536psin^' 

30473 660-C2 9697 cos 2-0 2581 cos 4-0 



4 

cos tjj 



1105920 



92160 



184320 



573 cos 6-0 121 cos 8-0 13 cos 10?/; (ir — 2-0) sin 2-0 



81920 



73728 



49152 



!)() 



3. Conditions at tp — ty/2 
The metric (|77|) should be symmetric about the plane ip = 7r/2. As a result we find that 

b 



10037T 4 

103680 



6 



(89) 



4- Verification of the first law 



We calculate the entropy of the black hole 



S = 



2G 5 



1 



/r 2 e 2 /37^ 4 „ . 

- + e Uo- + 3C2 



We compare the zeroth law T = k/(2tt) and the first law T" 1 = dS/dM to find that 

23tt 4 



C2 



1920 



To summarize, the location of the horizon, the entropy, and the temperature are 

7r 2 e 237r 4 e 2 " 



Ph = (Jt 
S 



1 



24 



tt 3 /* 3 r 



T = 



2G 5 
1 

277/1 



14 
1 - 



1920 

7r 2 e ^ 4 e 2 



8 384 
57r 2 e 437r 4 e 2 



24 



1152 



The near metric is given by 



B 



N 



P 2 H 



7r 2 e 



6 ^ COS ^ + 3]+ 5W 



[1570/i 2 + 17p 2 - 4(14p 2 - 395p, 2 ) cos 2-0 



+ (52p 2 + 580p, 2 ) cos 4-0 - 4(2p 2 - 25/i 2 ) cos 6^ - {p 2 - 2p 2 ) cos 8^] 



(90) 

(91) 

(92) 
(93) 
(94) 



(95a) 



,N 



jjN = 



l + p[(3p 2 ^ 2 )cos 4 ^l] + ^ 

+8(135p 4 + 88p 2 p 2 - 55p 4 ) cos4i/; + 4(45p 4 + 26p 2 p 2 - 15p 4 ) cos 6ip + 5^ 2 (4p 2 - /1 2 ) cos 8^] , (95b) 

7r 4 e 2 sin 2 -0 



[l800p 4 + 1324p 2 p 2 - 307/i 4 + 4(675p 4 + 422p 2 p 2 - 325p 4 ) cos 2^ 
+ 26p 2 p 2 - 15p 4 ) cos 6-0 + I 
[-540p 4 - 970/?p 2 + 385// 



V N = - 



l +1F ( P -p)sin tp cos ^ 144QAf4 

-5(144p 4 + 229p 2 p 2 - 385/i 4 )cos2V' - 5(p 2 - p 2 )(36p 2 + 82p 2 ) cos4^ - 55p 2 (p 2 - p 2 )cos6V>] , (95c) 
7r 2 ep 3 sin-0 cos 3 -0 7r 4 e 2 p 3 sin 2-0 



fi 2 2880/t 4 
2(90p 2 + 89p 2 ) cos 4t/> + 25/i 2 cos &ip] 



[990p 2 + 540p 2 + (720p 2 + 1207p 2 ) cos 2^ 



(95d) 
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VII. SUMMARY AND DISCUSSION 



In this paper we have calculated the metric of a small black hole in a five dimensional cylinder. The metric is found 



perturbatively, with an expansion parameter e = /i L = 8G§M(3irL 



M is the physical ADM mass and L is 



the circumference of the compact dimension, both measured at infinity. We calculate the metric, up to second order 
in e (fourth order in the ratio p/L). The asymptotic solution is given in Eas. (|12ll3ll7|) and in the appendix \K\ The 
near solution is provided in Eqs. (|21I95[) . 

We obtained the following expressions for the entropy and the temperature of the black hole 



S = 



T = 



2 L 3 



2G 5 
1 

2ixp 



3/2 



1 - 



1 



8 384 
57r 2 e 437r 4 e 2 ~ 



24 



1152 



(96) 
(97) 



These, and other results, presented below, were previously obtained in first order of e by Gorbonos and Kol 0] and 
by Harmark |17| . Using the Smarr formula we find the relative binding energy 



, . ir 2 e ttV 



(98) 



We have chosen the coordinate system such that the horizon is independent of the angle, ip, between the three 
dimensional space and the w axis. As the scale of the five dimensional radial coordinate, p, is fixed by the term of 
the metric, p 2 sin 2 ipdQ 2 , the radius of the horizon is uniquely determined. The horizon is located at p = pu as in 
Ed.ipHJl. We find that 



PH = P 



1 



237r 4 e 2 



24 



1920 



(99) 



A. Discussion of the black hole - black string transition 



The question of transition from black hole to black string, as a function of the mass parameter, has been discussed 
at length in the literature @ Q ■ One scenario proposed was the transition at the intersection of the black 
hole and nonuniform black string lines in the (M,n(M)) (mass-relative tension) phase diagram. We use the second 
order results for the small black hole to extrapolate to such a transition. We study three aspects of the transition; 
the mass-relative tension phase diagram, comparison of entropies, and the change of topology (using two different 
approaches) . 

A uniform black string is described by the metric 



ds 



BS 



2G 5 M 
Lr 



dt 2 



1 - 



2G n M 



Lr 



dr 2 



- 2 dQ 2 + dw 2 



The relative tension is constant, ubs = l/d—3 = 1/2, and the entropy is 



S 



BS 



4G 5 



-A-kL 



2G 5 M 
L 



3t3 



9vr J i 
~16G~5 



(100) 



(101) 



The entropy and relative tension of the non- uniform black string were have not been found numerically yet.jjj. Such 
a configuration exists for masses larger than the Gregory-Laflamme mass, which is, in terms of our variables, e > egl- 

In Fig.@ we draw the (e,ro) phase diagram for the uniform black string, the non-uniform black string (just the 
leading order term near the GL transition), and the black hole branches. The black hole turns into a nonuniform black 
string beyond the transition point, e = 0.1. We did not include recent numerical results 0, which, due to numerical 
difficulties, have large errors |2l| 

In Fig.QJ we draw the entropies of the black objects in units of L 3 /G$. Comparing (|96|l and l|101|) we find that 
the enropies of the uniform black string and the black hole are equal for cbh-bs — 0.102. Below that value the black 
hole has higher entropy, and above that value the black string has higher entropy. This simple check suggests that the 
black hole will be unstable for e > 0.102. Note again that our second order expansion in e shows good convergence. 
Using the first order expansion only the intersection of entropies would occur at e = 1.01. 
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FIG. 3: Phase diagram for the uniform black string (dashed green line), non-uniform black string (blue line), and black hole 
branches (solid red continued in dashed red in the region where the black hole turns into a nonuniform black string). 
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FIG. 4: The entropy (in units of L A /G5) of the uniform black string (dashed green line), non-uniform black string ( dashed red 
line), and black hole branches (solid red line). 



A transition between a black hole and a black string requires a topology change of the horizon. This change is likely 
to happen when the horizon of the black hole fills the compact dimension and touch itself. At this point a small pinch 
of the horizon will turn it into a non- uniform black string. In the notation we use, Ea. Hll[> . this will happen when 

L , , 

p H = - e„ = 0.0957 . (102) 

The first order formula would give e p = 0.094. Again, we see a fairly good convergence, making our result reasonably 
reliable. We will see a further indication for this point in another, independent, calculation of the critical value of e. 

The apparent agreement of e p with £bh-bs (equality of the entropies of the uniform black string and the black 
hole) is somewhat problematic. We see no mathematical reason for such a agreement, though it would be interesting 
if these points also coincided. Note that for geometric reasons the black hole branch should cease to exist above the 
critical value of e so the stable state of the system is the uniform black string, which has a higher entropy than the 
nonuniform black string. 
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B. Discussion of the shape of the horizon 



In terms of the coordinates we chose the horizon is located at constant p. However, this does not mean it is spherical. 
Compactification breaks the 0(4) symmetry. The horizon can be an oblate (elongated perpendicular to the compact 
direction) or a prolate (elongated parallel to the compact dimension) ellipsoid. There is no generally accepted method 
for distinguishing between the oblate and the prolate configurations. Some authors [1,^2 define the 'eccentricity' of 
the black hole as 

e = 4l - 1 , (103) 

where A\\ is the maximal area of a cross section of the horizon parallel to the compact dimension, and A± is the 
maximal area of a cross section of the horizon perpendicular to the compact dimension. 



A,, = J 7r dip J" dtj> p 2 sin tp sin 6^U(p, if>) =^Ph / # sin ipy/U(p H ,ip) , (104) 

P=Ph ,S=7r/2 







A x = d0 J dcj> p 2 sin 2 ip sin 6 = 4irp 2 H . (105) 

p=PH ,l/>=7r/2 

A prolate (oblate) horizon has positive (negative) eccentricity. We use the near solution (|77|l to find that e = 47r 4 e 2 /135. 
We see that the horizon is prolate. 

Another possible measure of the eccentricity is given by the intrinsic Ricci scalar of the horizon, which is calculated 
in terms of the three metric 

ds 2 H = p 2 H [U{p H ,ij)dil) 2 +sin 2 V (d6 2 + sin 2 6>d</> 2 )] . (106) 
Using the near solution, the intrinsic curvature of l|106|) is 



R (3) = 4- 

Ph 



27r 4 e 4 , 

-^-(1 + 5 cos 2V0 



(107) 




It is maximal at ip = 0, which indicates that the horizon is indeed prolate. 

One might speculate that a prolate horizon will tend to grow in the periodic direction and turn into a nonuniform 
black string. However, the circumference of the compact direction grows as well 12] and the transition into a black 
string is not obvious. In it is suggested to measure the 'inter polar distance', which is the proper distance between 
the poles of the prolate horizon along the compact dimension. In principle, to calculate this distance we need to break 
the integral defining this distance into two parts. In the first part we should use the near solution and in the other 
we should use the asymptotic solution. In the near horizon region one should integrate along ip = from p — pn to 
some arbitrary point p — z. In the asymptotic region we should integrate along r — from w = w(z) to w = L/2 

dw . (108) 

arcsin irz/ L 

In |l2j | it was found that the zeroth order approximation is L po i es = L, which means that the circumference grows 
enough to make room for the black hole. 

If the mass (e) value, at which the black hole fills the compact dimension, is small enough then it is possible that 
the integral over the near solution is sufficient in (|108f) . or, in other words, we can choose z = L. Then we obtain 

^^149^ h003^ 

L 60tt 24 1440 V ; V ' 

First of all, we see that the distance between the poles of the horizon is a decreasing function of e. If we solve the 
equation L po ios = we obtain e cc = 0.146, which is in a rough agreement with values we obtained for the critical 
value, e c , above. Taking the first order term only we would obtain e cc = 0.355. Clearly, the convergence of e cc is 
not as good as those of e c and e p . One reason for the poorer convergence is that for very small e (|109fl must fail, as 
a significant contribution coming from the asymptotic solution was omitted. It should work, however in the region 
where L po i cs <C L. Fortunately, in this region the e-expansion is still fairly reliable and the rough agreement of e c ~ e cc 
is encouraging. 
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APPENDIX A: SECOND ORDER ASYMPTOTIC SOLUTION 



The second order asymptotic solution is given by Eas. l|t)9|l 



i {2) 

Hi 



1 {H) + h {2) 

l tt T lift , 



2'Ht,r ' n rr 1 
ll ww — 2 U' r ' w WW ' 



h (2) = W (2) + ~ h (2) 
rw ,r 1 rw 

The the components of the inhomogeneity, h^ 2 \ are 



;( 2 ) - 



~h 2 + Wh w , 



~ h ?r = T^[(h,r) 2 + (h,w) 2 ]-^[hh, r + 4Wh, rw } + (W, r ) 2 -rK^ 



U 2) 
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T — [(h r ) 2 + (h w ) 2 ] - ~ [2W, w h r + W h <rw ] + -^h 2 + (W, w ) 2 -rK,, 

K-^[hh :W + 8W, r h,r] + W !W W, r , 



where h = h[]\ W = W^ l \ and K, ww — hh^ rw /8 — h tW h jr /4, 

h = w 



ttL sinh ^ 



r(cosh ^ - cos ' 



W = W 
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7rLsin22H 



irL 3 2nr 
K = 16^ COth — 



12(cosh^ -cos^iH) 



2nr 2nr 
1 — coth — — - 



3L . 2itr 2-kw 

smh 2 cos h 5 

2-kt L L 

( irw , irr\ irL 3 sin M 22Z coth 
arctan tan — coth — + L L ' - 



2jrr 



2-nw 1 
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(Ala) 
(Alb) 

(Ale) 
(Aid) 

(A2a) 

(A2b) 

(A2c) 
(A2d) 

(A3) 
(A4) 
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The function was chosen such that W(r ^S> L,w) — ©(r 1 ). It satisfies Eqs. (|41I49|I . and it has a structure 



similar to h[]\ 

(2) 

K is an odd function of w, but K (r, w = L/2) ^ 0. However, ^-symmetry implies that h r w{r, w = L/2) = 0. This 
means that the gauge function, W^ 2 \ has nontrivial boundary conditions at w = ±L/2 



WW(r,w = L/2) = -WW(r,w = -L/2) 



ir 2 L 3 , 2irr 

=- coth 

32r 3 L 



2nr 2irr 

1 coth 

L L 



(A6) 



To match the asymptotic and the near solution we have to transform the asymptotic solution to the (p, ip) coordi- 
nates, using (|llf> . and expand to order L~ 4 . We also expand the gauge function as 



W^ 2 \p, ip) = L 4 



(A7) 
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The asymptotic metric in the (p, 0) coordinates, expanded to order L 4 is 
,2 n 2 /! 2 sin 2 ip n 4 p 2 p 2 sin 4 
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+ 
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(A8a) 



(A8b) 
(A8c) 



L 2 



[p 3 sin cos 3 + ^ 2 psin0cos 3 0(sin 2 + 2) + p 2 psimpW 2 , p (p, 0) - p 2 cosipW 2 ^(p, 0)] . (A8d) 
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